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ABSTRACT 



We study an ideal Bose gas of N atoms contained in a box formed by two 
identical planar and parallel surfaces S , enclosed by a mantle of height 
a perpendicular to them. Calling r the mean atomic distance, we assume 
S ^> Tq , while a may be comparable to r^. In the bidimensional limit 
(a/r = 7 << 1) we find a macroscopic number of atoms in the condensate at 
temperatures T ~ 0(1/ log(iV)); therefore, condensation cannot be described 
in terms of intensive quantities; in addition, it occurs at temperatures not too 
low in comparison to the tridimensional case. When condensation is present 
we also find a macroscopic occupation of the low-lying excited states. In 
addition, the condensation phenomenon is sensitive to the shape of S. The 
former two effects are significant for a nanoscopic system. The tridimensional 
limit is slowly attained for increasing 7, roughly at 7 ~ 10 2 — 10 3 . 
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1 INTRODUCTION 

A system of N non-interacting particles is in a condensed phase if a non- 
negligible fraction of them lie in the ground state; however, for small values 
of N , say 10 3 10 5 , the condensation boundary becomes blurred. As it 
is well-known, for an ideal gas obeying Boltzmann statistics, condensation 
occurs at vanishingly small temperatures, irrespective of the dimensional- 
ity d of the system ( T ~ 0{1/N 2 } for d = 1; T ~ 0{1/N} for d = 2 
and T ~ 0{1/N 2 ^ 3 } for d = 3). In contrast, for Bose statistics, conden- 
sation begins at a finite temperature in three-dimensional systems, while it 
is absent in one- and two-dimensional systems at least in the mathema- 
tial limit N -> oo ( T ~ 0{l/N} for d = 1 and T ~ £>{!/ log(iV)} for 



d = 2, see Ref. fllQfl). The presence of an strong enough external potential 
can, however, bring about a Bose-Einstein condensation (BEC) in two and 
even one dimensions, at a finite temperature 0. 

Since Bose and Einstein predicted the condensation phenomenon, there 
has been a great interest into achieving such a state in an actual system. The 
existence of superfluidity in 4 He has been considered a sort of corroboration 
of BEC. However, superfluid Helium is a poor example of a non-interacting 
particle system, since the mean atomic separation is comparable with the range 
of the atomic potential. Therefore, a complete explanation of superfluidity 
must rely in Many-Body Theory. 

The goal of obtaining BEC in a gas of nearly non-interacting atoms was 
recently achieved by Anderson et al. who combined the techniques of laser 
and evaporative cooling of magnetically trapped atoms, implementing the lat- 
ter one by means of an ingenious "time orbiting potential" magnetic trap. 
They obtained temperatures of ~ 10~ 7 K for a sample of a few thousand 
atoms trapped in a region of linear dimensions ~ m and were thus able 
to achieve BEC. The nanoscopic character of the experiment motivates a more 
detailed study of BEC for such small systems. 

The purpose of this work is to analyze the condensation process of an ideal 
Bose £f function of the dimensionality of the system and the shape of 

the container. Since a rigorous infinite two-dimensional gas is not a realistic 
system, we shall consider here a true tridimensional gas enclosed in a finite 
container; our aim is to understand how the properties of the gas change 
as the container is shrunken in one of its dimensions (say, z), expecting a 
continuous passage from tri to bidimensional behavior. In fact, at enough small 
temperatures, when the de Broglie wavelenght A(T) becomes smaller than the 
container width, there is not enough thermal energy to excite the degree of 
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freedom associated to the z direction. Therefore, one dimension "freezes" 
and our system behaves like an "effective" two-dimensional gas, in spite of 
its actual tridimensional nature. If the container width is comparable to the 
interatomic mean distance, condensation must occur after the bidimensional 
regime has been established. 

Interest in this problem dates back to London's suggestion to use BEC as 
a qualitative model for understanding superfluidity in He 4 [|J. Experiments 
showed that He 4 remained a superfluid when forced to flow inside narrow 
capillaries. On the other hand, superconductivity was shown to persist in thin 
films. This prompted the study of "superproperties" in finite geometries. In 
one of the first studies @, an ideal bose gas was enclosed in a box of finite 
volume (L x D x D) and the number of particles in the ground state was 
computed numerically as a function of the temperature. Large deviations from 
the bulk result were found in the one-dimensional limit (L » D). Attempts 
to ascertain the dependence of the superfluid transition temperature for He 4 
on the dimensions of narrow channels §, were followed by discussions on the 
ambiguity of the concept of a transition temperature in a finite geometry, 
where a rigorous phase transition is absent [|7|, §. Vanishing of quasiaverages 
in partially finite geometries (one or more dimensions finite while one or more 
dimensions extend to infinity) led to the suggestion IJ that these geometries 
were not good approximations to thin films and pores found in laboratories. 
For a strictly finite geometry, BEC was shown to occur for an ideal gas||. 

While in most of these earlier works, the goal was to achieve a qualitative 
understanding of superfluid properties, we are now genuinely concerned with 
describing the behavior of a finite number of nearly-independent bosons, in- 
side a confining potential, such as in the recent experiments mentioned at the 
beginning0. The recent advances in nanotechnology and laser optics allow 
researchers now to manipulate molecules and their environment at the meso- 
scopic and nanoscopic levels, rendering the investigation of quantum effects in 
a system with a finite number of components (i.e., away from the thermody- 
namic limit), an interesting field of study. 

In the present work we shall address the problem of crossover between two 
and tridimensional behaviour in BEC. We shall give both, exact numerical and 
asymptotic analytical results. For the sake of completeness we shall refresh 
some old results. However, our anlytical expressions will give a compact and 
clear sight of the phenomenology under consideration; in particular, they show 
the non-intensive character of BEC near the two-dimensional limit, giving 
simple estimations for the temperature of condensation (we remember that 
the latter temperature is not sharply defined in this case). We shall also 
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clarify the effect of geometry in BEC, which can be summarized by some 
characteristic constants associated with the shape of the container, listing them 
for some particular cases; incidentally, the different boundary conditions will 
be considered like peculiarities of the topology of the container, affecting the 
values of former constants. The role of lower excited states and its macroscopic 
occupance shall be also addressed. Finally, we shall also analyze the continuous 
transition between two and tridimensional geometries. 

2 THE MODEL 

We study an ideal Bose gas of iV atoms contained in a box formed by two 
identical planar and parallel surfaces S (say, in the x, y plane), enclosed by 
a mantle of height a perpendicular to them. We call ro the mean atomic 
distance, and introduce the parameter 7 = a/ro . Since the container vol- 
ume is Sa , it holds 7q = Sa/N . We assume that S ^> a 2 (or equivalently 
N ^> 7 3 ), while a may be comparable to ro- 

We introduce a characteristic temperature for the ideal gas, associated to 
the energy of a particle whose de Broglie wavelenght A(T) is comparable to 
the mean intermolecular distance r , 

h 2 

^° 2irmkBro 

where m is the mass of a particle and h , kg are the Planck and Boltz- 
mann constants, respectively. For T < Tq the bidimensional behaviour can 
be associated with 7 ~ 0(1) , while 7 0(1) corresponds to the tridimen- 
sional limit. We also define the "reduced" temperature r = T/T ; now the 
(tridimensional) BEC critical temperature is r c = 0.527201068760... ; see 
Ref. 0. 

The energy dispersion relation for a particle enclosed there can be expressed 

like 

E(v,£) = k B T{ V C + a£ 2 } (2) 

here o = 7r/[r7 2 (g + l) 2 ] , where q = for ciclic boundary conditions (CBC) 
and q — 1 for impenetrable boundary conditions (IBC) in the z direc- 
tion respectively; £ = 0, ±1, ±2, . . . for CBC, while £ = 1, 2, 3 ... for IBC. 
7] = tt/tN s , where N s = N/j = S/r 2 , is the number of particles in a neigh- 
bourhood r of container base 5*. The numbers j^} describe the spectrum 
of the bidimensional Helmholtz equation for the surface S; they depend only 
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on the shape of S and the boundary conditions used for the x, y plane. We 
choose the ground state as energy zero; therefore we shift ^ — > Q — £q , while 
£ 2 — > £ 2 — 1 for IBC in the z direction. 

We shall use the asymptotic condition for the bidimensional energy 

D(0 d£ = {Number of £ v e [£ ,£ + d£\} = (2< - C s ) d£ (3) 

valid for £ d£ 3> 1 . Here C s is a constant that depends of the shape 
of S ; therefore, the contribution — C s in Eq. ([!]) represents the finite size 
corrections to the density of states [llj ; for example, C s = for a square with 



CBC, C s = 2 for a square with IBC, and C s = y/n for a circle with IBC. 
The mean occupation of the state E(i/,£) is given by 

A/>,£) = N-n(u,£) = {Cexpf^ + a (f - q)] - . (4) 

here ( = exp(— /i/A^T) , where /i is the chemical potential. Ac- 
cording to the former relation, the ground state mean occupation is 

Af = N-n = !/((- I) ■ 

For IBC in the z direction we shift the index £ — > £ — 1 , in order to 
associate £ = to the ground state. Due to the fact that iV ^> 7 3 , it holds 
that n(i/, £ 0) ~ 0{1/N}, while, for £ = the lower excitations in the 

{x, y} plane may have a macroscopic population. In fact, when a macroscopic 
amount of particles lie in the ground state, relation (^) leads to 

n&,0) « (- + — (^) 2 + ^(^)) 1 (5) 
I no r TV J 

The equation for ( is 1 = Y, v i n(u,£). In order to study the conden- 
sation phenomenon, it is necessary to separate the former sum into contri- 
butions from states with macroscopic ( n(u,£ = 0) ~ 0(1) ) and microscopic 
( n(u,£ ^ 0) ~ 0(1/N) ) occupancy. Only in the latter contributions we can 
safely approximate J2v ~^ I d£D(£) since then n(u,£) is a smooth function 
of v. But the macroscopic contributions n(i/, 0) must be summed explicitly, 
at least for the lower energy states. Over a high enough value of v it is safe 
to replace this sum by an integral; defining < R < , R ^> 1 , the 
equation for ( becomes 



1 Y) + 5>&"°) 



*(c .... 

+ {q + 1)N J2j d£D(0{(exp[ V e + o-£(£ + 2q)} - l}~ + 



7 



- j R de^(e)[cex P (^ 2 



-i 



(6) 



where 1/[N {Q — 1)] = n is the fractional ground state occupation. The 
bulk contribution to D(£) can be integrated in a closed form; the finite size 
correction to D(£) gives a negligible contribution, except for the case when 
condensation has begun and £ = ; in the latter case we use Eq. fl5|). Thus, 
the equation for ( takes the form 



1 



N(( - 1 ) 
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+ ^E[Cexp(^ 2 )-l 
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(7) 



This expression is suitable for numerical computation. In practice, we 
have summed over the L = 9000 lower states in the case of a rectangular 
base, and L = 1241 states for the circle; for such values of R the error on 
approximating the sum by an integral is negligible. In addition, for a given 
L , the specific choice of R in the range £l < R < is a ls° irrelevant 

(the associated variation in n is lower than 0.4 x 10~ 4 ). 

Some Approximate Results 

Before showing and discussing the numerical results, it is useful to obtain 
some asymptotic relations for the case no 3> 1/N and a not too large value 
for 7 , say, a = 7t/[t7 2 (1 + q) 2 } > l/[jVra ]. Then, it follows from Eq. (g), 
( = 1 + l/(Nn Q ) 1 . In addition, replacing Eq. (|5|) into Eq. (|7|) and taking 
the limit R — ► oo we obtain 
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where 



c\ oo 

- ^log{l - exp[- r£(£ + 2g)]} (8) 
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(i) Small Temperature Case— We shall assume that condensation is 
well established, and that the inequality n ^> r/7 holds. We expand in a 
Taylor series the second term of Eq. @ , concluding the following relation for 

n 
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Solving for n we obtain 
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where m(r) = 1 [log(arA?" s ) + $(cx) ] (10) 

7 

An estimation of the accuracy of this approximation is 
An = n — {no} exact ~ CV 3 /(7 3 nQ) , where C ~ 0.2 , at least for square 
and circle geometries. Thus, the error strongly decreases as n , N s or 7 
increases. However, the error is very small for n = 0.2, where we have 
An < 0.01 if > 10 4 , 7 > 1. When N > 10 s the error An < 0.003 , 
even for a very small amount of condensate, uq > 0.08 . The accuracy of 
Eq. (|i0|) decreases if the shape of S is very "oblong". It is important to 
remark that the constants a and k only depend on the shape of surface 
S ; we include a table with the values of these constants for a circle and a 
rectangle b\ x 6 2 ; in the latter case we use different boundary conditions and 
values for the ratio p = b\jb% (e.g. p = 1 corresponds to a square). 

As a consequence of these geometric constants, we conclude that condensa- 
tion is inhibited when the shape of S is very oblong; the latter one is consistent 
with the absence of condensation in a nearly one-dimensional container Jl| , in 
accordance with our exact calculations (see Fig. 1). 

(ii) Intermediate Temperature Case— Using Eq. (||) in the case 
1/iV < % < 1 , when condensation is scarcely incipient, we have concluded 
the asymptotic relation 
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exp 



-27(1 + q)\- 

7T 



(11) 



where 77 = 7/ log(iV) roughly corresponds to the temperature threshold for 
condensation appearance. In this way, the typical temperature for condensa- 
tion in a quasi-bidimensional system roughly corresponds to tj ~ r c / \og(N) , 
where r c is the tridimensional critical temperature for BEC; the latter one 
implies that, typicaly, the quasi-bidimensional condensation occurs at tem- 
peratures one order of magnitude below the tridimensional case. 

Results and Discussion 

We have solved Eq. (0), determining uq{t) for several values of N , 7 and 
different shapes of S: a circle, rectangles with different flatness p . We have 
always used IBC in the normal direction with respect to S. We display our 
numerical results in figs. 1-4. In what follows we summarize the principal 
consequences of these results. 



(i) The shape of the base S influences the condensation phenomenon, 
speciality in the nearly bidimentional limit, 7 ~ O {1} . This effect is also 
enhanced for relatively small values of iV . In Fig. 1 we illustrate this effect 
for 7 = 1 and N = 10000 . We see that a IBC-square base yields the larger 
values of n (r) (at least for not too small values of no); however, the circle 
gives only slightly lower values for n in comparison to the IBC-square, while 
the CBC-square curve Uq(t) lies below the two former ones. Finally the lower 
values of tiq{t) in Fig. 1 correspond to IBC-rectangles; these curves descend 
on increasing the flatness p , as the system approaches to a one-dimensional 
behaviour. The latter one is in accordance with the absence of BEC in that 
limit §. 

This "shape effect" on no(r) is important in the 2-d limit, diminishing 
upon increasing 7 or N. Therefore, it can only be observed in truly nanoscopic 
systems. In the following item we shall explain this behaviour. 

(ii) According to Eq. the low-lying excited states also have a 
macroscopic occupancy once condensation begins. In Fig. 2 we display the 
occupancy of the lower seven excited levels in the case of circular base, 

N = 10000 and 7 = 1. The associated quantum numbers are v = (m,j) = 
(0,1), (1,1), (2,1), (3,1), (1,2), (0,2) and (2,2). Here m is associated to the 
axial component of angular momentum L z , while j enumerates the eigenen- 
ergies for a given m. Obviously, while m = is not degenerate, m^fl has 
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degenerancy two. The effect of degeneracy becomes evident in the asymptotic 
behaviour r — > oo , where the occupancy of non-degenerate states is one-half 
of the degenerate ones. The maximal occupancy of an excited level corresponds 
to (1,1) (say, the first excited state), giving nfi^TMax — 0.13) = 0.093 . 

In the general case, the macroscopic occupancy of excited levels is enhanced 
as 7 or JV degreases. An approximate relation for the maximal occupation 
of these excited levels is given by 



n u (r Max ) ps ^ with 9 = n log 



Na 



(12) 



_log (N s a) 

and r Max ~ ; here D u is the degeneracy. Thus, this maximal occu- 

pancy depends on N like l/log(iV) (thus increasing as N decreases, being 
appreciable in truly nanoscopic systems), and it is inversely proportional to 
the bidimensional excitation energy E(y, 0) — E(0,0) (see Eq. (0) ). We re- 
mark that Eq. (|12"D is valid under the condition 7 ~ 1 ; for larger values of 
7 the occupancy of the excited levels strongly decreases (for example, us- 
ing the same N s of Fig. 2, but 7 = 10 , we obtain a maximal occupancy 
n (1 ,i)(T Max = 0.47) = 0.037 ). 

According to relation fll2|) , the macroscopic occupancy of the lower excited 
levels is strongly affected by the particular sequence , which, in turns, 

depends of the shape of S . The latter one implies that the occupancy of the 
ground state is also affected by the base shape. In particular, the amount of 
condensate no increases if the gap between the ground state and the first 
excited levels increases, i.e., when the level density decreases. For a fixed base 
area but different shapes, we can modify the former magnitudes and therefore 
modify n . 

In order to obtain a further feeling for the effect of the shape in the amount 
of condensate no(r) , it is very illustrative to note the following numerical 
result: The sum of occupancies of the ground state and the first excited level 
is roughly the same for the circle, the CBC-square and the IBC-square; the 
latter one, near the theshold temperature of condensation. Therefore, the 
shape effect on uq roughly corresponds to the difference in the populations of 
the first excited levels. Although this estimation is very crude, and certainly 
not applicable to arbitrary shapes (specially if one of them yield a quasi- 
degenerate first excited level), it gives a first estimate for the shape effect on 
the thermodynamics of the condensate. 

The effect of shape in condensation is also evident from relations (^) 
and (|T0|), since the geometrical constant a increases as the density of levels 
in the neighbourhood of ground state decreases (see Eq. (|10D ). An increase 
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in a implies a decrease of n {r) , in accordance with Eq. (|10D. 

The results of Fig. 1 is consistent with this explanation. In addition, since 
the occupancy of the low lying excited levels is only important for small values 
of 7 or N, the effect of shape in the ground state occupanncy uq(t) is ap- 
preciable under the same conditions. Therefore, this "shape effect" can only 
be observed in truly nanoscopic systems. 

The latter observation is consistent with the relation between the density 
of states D(E) an the system dimensionality. For example, in a fully tridi- 
mensional system, D(E) ~ y/~E , vanishing as E — > ; this low density is 
responsible of a genuine BEC at a relatively hight temperature. In the bidi- 
mensional case, D(E) ~ {Constant} ; while the onedimensional case corre- 
sponds to D(E) ~ l/y/E , diverging as E — > . Accordingly, the interme- 
diate density of states of the bidimensional case is consistent with condensation 
phenomenon at temperatures one order of magnitude below the tridimensional 
case. Finally, in the one-dimensional case, the high density of states near the 
ground state is in accordance with the absence of condensation in that case. 

(Hi) According to relation (|II]), the condensation phenomenon is highly 
non-intensive in the case of a small 7 , due to the dependence of Uq{t) 
on log(iV) . Moreover, using that relation we can estimate an "effective crit- 
ical temperature" r e ff ~ 7/ log [a AT/ log(iV)] ~ 7/log(iV) ; this estimation 
is consistent with 77 of Eq. (|TT|) and Tuax of Eq. (fl"2|). 

In figs. 3 and 4 we plot no(r) for iV = 10 4 and 10 8 for the case of a circle; 
in each figure we take different values of 7. In general terms, these figures 
confirm that condensation is inhibited on decreasing 7 or increasing N. 
Inspection of these figures confirm the non-intensive character of n (r) . For 
example, the curves 7 = 1 of Fig. 3 and 7 = 2 of Fig. 4 nearly coincide, in 
agreement with the fact that r e ff ~ j/log(N) has the same value in both 
cases. 

The condensation phenomenon has a nearly 2-d character for small values 
of 7 , as long as the excitations in the ( "thin" ) z direction are very improbable 
at temperatures compatible with condensation of a macroscopic number of 
particles. An upper bound for this 2-d character can be roughly estimated as 

7 < 2 for a nanoscopic system (N s ~ 10 4 <-> 10 s ). The tridimensional limit 
is slowly attained as 7 increases; roughly, for 7 ~ 10 if N = 10 4 , and 

7 ~ 50 if N = 10 8 . In this way, on increasing 7 , the tridimensional limit is 
attained faster for a small number of particles; this result seems natural, since 
for a fixed value of 7 the base width (~ y/S) and container height (a) may 
be comparable if N is not too large, but (y/S ^> a) for a greater value of N. 
Another important results illustrated in Fig. 3 are: (a) the finite size 
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effects, which yield a larger amount of condensate for 7 = 20 in comparison to 
the tridimensional case, (b) The condensation phenomenon begins gradually 
for low values of N and 7 , in close accordance with the qualitative result (|TT|); 
while larger values of N and 7 yield a more sharp phase transition. 
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Figure Captions 

Figure 1: Number of particles in the ground state as a function of the reduced 
temperature for several different shapes of the container base. In all cases, we 
have IBC, except for the lower square, where CBC are used. (7 = 1, N — 10 4 ). 

Figure 2: Occupancy of the lower seven excited states for the case of a circular 
base. Each state is labelled by (m,j) where m is associated to L z while j 
enumerates the eigenenergies corresponding to a given m. (7 = 1, N — 10 4 ). 

Figure 3: Circular base: Number of particles in the ground state as a function 
of the reduced temperature, for different values of 7. (N = 10 4 ). 

Figure 4: Same as in fig.3, but for N = 10 8 . 
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